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In a framework for long-range density-functional theory 
we present a unified full-field treatment of the asymptotic 
van der Waals interaction for atoms, molecules, surfaces, and 
other objects. The only input needed consists of the electron 
densities of the interacting fragments and the static polariz- 
ability or the static image plane, which can be easily evalu- 
ated in a ground-state density-functional calculation for each 
fragment. Results for separated atoms, molecules, and for 
atoms/molecules outside surfaces are in agreement with those 
of other, more elaborate, calculations. 

PACS numbers: 71.15.Mb,31.15.Ew,34.20.-b,34.50.Dy 



The ubiquitous van der Waals interaction needs an ef- 
ficient and accurate description in many contexts such 
as interacting noble-gas atoms, van der Waals com- 
plexes, physisorption, interacting macroscopic neutral 
bodies, liquid-crystal interactions, solute-solvent inter- 
actions, and soft-condensed matter. Epr dense mat- 
ter the density Jmictional theory (DFT),lil with its local- 
density (LDiUB^ and generalized-gradient approxima- 
tions (GGA)j3i3 is a clear success. Ground-state and 
thermodynamic properties of increasingly more complex 
systems are now being calculated with a practically very 
useful accuracy. As the world around contains far more 
objects than just hard solids, a generalization of these 
methods to also account for the van der Waals forces is 
in great demand. -These forces are an inherent property 
of the exact DFTQ and it is thus a question of providing 
an approximate van der Waals density functional that 
is generally applicable, efficient, and accurate, and that 
pep. definition is ar-fiiHctional of the density only. Earlier 
wea't3 and othersU^t3 have proposed such functionals and 
shown them to give useful results in significant applica- 
tions. For a review, see Ref. Until now, however, 
there has been a certain asymmetry in our treatment of 
"small" and "large" objects, respectively, which will be 
remedied in this paper. 

The starting point for our functionals is the exact ex- 
pression for the exchange-correlation energy _Exc as an 
integral over the coupling constanlj^^, the so-called adi- 
abatic connection formula (ACF)tlo 

E^Jn] = - f d?rd\' X 
^ ^ 2 / r-r' 



where h = n — n, n being the density operator, and 
(. . ■)n,\ means that the integration is performed with a 
potential V\ present, keeping the density equal to n(r). 
To second order perturbation theory in the interaction 
Vab between two separated objects a and b the ACF can 
be cast into the formoEj 



V;fc(R -I- ri - r2)Kh(R + rs - 1-4) x 

/ — na(ri,r3;iu)nfc(r2,r4;m) . (2) 
Jo 

Our evaluation of this energy is based on two approxima- 
tions. First we introduce a local dielectric function that 
depends on the local electron density. Second, we limit 
the volumes of the interacting objects by using a cutO|fiE- 
an idea first introduced by Rapcewicz and Ashcroft,E£l 
outside which the response to an electric field is defined 
to be zero. We thus have the dielectric function 



e((jj;n(r)) = 1 — K{n(r))- 



i(n(r)) 



where 



Wp(n(r)) = 47re^n(r)/m . 



(3) 



(4) 



The cutoff is implemented via the function K(n(r)), which 
is either unity or zero, following the notion discussed ear- 
lierO that the local approximation for dielectric response 
greatly exaggerates the response in the low-density tails, 
where it is better to assume no response at all. 

These approximations are common for all of ouftrtiiree 
model systems - interacting atoms or molecules, BEJLif*i| 
atom or molecule interacting with ajplanar surface |lJ'E£l 
and finally two interacting surfaces.llB However, we have 
earlier treated the electrodynamics on different levels for 
"small" and "large" objects. Normally, local electrody- 
namics means a local relationship between the polariza- 
tion P and the total electric field E 



P(r,c.)-^[e(^;n(r))-l]E(r,c.), 



(5) 



dA [(n(r)n(r'))„,A-<5(r-r')Mr)}] , (1) 



which we use for surfaces. For atoms, however, the calcu- 
lations get somewhat more complicated than in the sur- 
face case when using Eq. (H). For instance, the electro- 
dynamics must be solved numerically for each frequency. 
Our earlier approach for atoms and molecules has ap- 
proximated the local polarization by 
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P(r,a.) = f [6(c.;n(r)) - i] EappUed(r a.) ^ 
47r e(a;; ?i(r)j 

which imphes the relation 

E(r, uj) = Eappiicd(r, '^)/e(t^; ri(r)) , (7) 

which is certainly wrong for macroscopic objects, 
but givesm surprisingly good results for atoms and 
molecules B'tZl However, the optimum cutoff position de- 
fined by K(n(r)) is somewhat larger for the approximate 
electrodynamics, and this fact seems to mitigate some of 
the deficiencies of this previous approximation. To ob- 
tain a unified treatment for different objects, and also to 
test our approximation for the dielectric function, this 
paper presents the electrodynamics using Eq. also for 
atoms and molecules. We apply it to the asymptotic van 
der Waals interaction of separated atoms, molecules, and 
parallel surfaces ,E3 and show the results to be in agree- 
ment with those of other, more elaborate, calculations. 

For two widely separated atoms or molecules the van 
der Waals energy is given by ~ —Ca/R^, where 

the van der Waals coefficient iaHj'EiJ 

3 f°° 

Cg = — / duai{iu)a2{iu) (8) 
Jo 

and aj{iu) is the polarizability of atom j. To calculate 
a{iu) we solve V • D(r) = 0, for each frequency in the 
presence of an applied constant electric field Eq. The 
displacement D(r) is given by e(ijj; n(r))E(r) and E(r) — 
— V0(r). Thus we solve the equation 

V-(e(w;n(r))V0(r)) =0, (9) 

where e{uj;n{r)) is given by Eqs. (^) and (^, and with 
the boundary condition that E approaches Eq at infinity. 

The cutoff function K(n(r)) must however first be de- 
fined. For surfaces the cutoff isJpund by taking the static 
limit of the surface response,E£l which implies that the 
cutoff should be defined by the static image plane d{Q). 
In order to have a common cutoff scheme for both atoms, 
molecules, and surfaces, and in addition to implement 
the requirement introduced for surfaces that the static 
polarization response be accurate, it is expedient to sim- 
plify the scheme used previously for atoms and molecules. 
This is done by choosing the cutoff function k according 
to 

K{r) ^ 0{n{r) - c) , (10) 

where c is a constant. 

For atoms and molecules, compared with the original 
scheme that uses both the density and its gradient, the 
practical effect of (|lO|) is to eliminate any cutoff in the 
intra-shell regions. We have found that inclusion of the 
cutoff in the intra-shell regions as before results in a me- 
dian reduction of the predicted values of the Cq coeffi- 
cients of 12% for the atom pairs calculated here. The 
extent to which these intra-shell corrections should be 



included even in principle is arguable, and since they are 
small, we have therefore opted for the simpler scheme 
(|To|). Adapting the analogue of the proceduxe used when 
the full-field scheme is applied to surfaces ,E£l we fix the 
constant c in ( p^ ) so that the static polarizibilities are ac- 
curate. For a spherically symmetric species, this means 
that the volume V inside which the step function in ( p^ ) 
is nonvanishing is simply V = (47r/3)a(0). For species 
without spherical symmetry, we choose c so that the 
isotropic polarizabilities, a(0) [see Eq. (|l^)], are correct. 
This scheme seems to underestimate the anisotropy of 
the molecular polarizability; if an accurate anisotropy is 
important the cutoff function should be modified so that 
the elements of the diagonalized static polarizability ten- 
sor are reproduced. 

The solution of Eq. (H) is dcjie with a finite element 
method with an adaptive net.E3 In this way we have a 
general method for all geometries. To secure a reasonable 
numerical accuracy (5%) at small frequencies, we here 
represent a(iu) by the expression a + b/ (1-|-m^/c^), where 
a, h and c are fixed by a smooth continuation of high- 
frequency results. In Table | calculated van der Waals 
coefficients for a number of pairs of identical atoms are 
given, together with the static polarizability used when 
defining the cutoff. In Figure |] calculated Cq values both 
for identical and mixed pairs of atoms are plotted against 
results from more accurate calculations. The values com- 
pare very well with results from other calculations, with 
the close agreement indicated by the narrow spread of 
the points around the diagonal. Especially the results for 
alkali and alkaline earth atoms are rtmch improved com- 
pared with our earlier calculations.B In Figures |^ and 
^ the dynamic polarizability from our calculations are 
compared with reference calculations. Figure || for He 
is a worst-case example, with a 12% error in Cg, while 
Figure ^ for Be is a best-case example, with a Cg right 
on the reference value. 

Results for a few molecules are given in Table ^ with 
results agreeing very well with literature values. The 
largest molecule for which we have so far calculated the 
van der Waals coefficient is fuUerene, Cgo. Recently the 
dispersion energy between two fuUerenes has been com- 
puted from first principles in time-dependent DFT, whi^ 
gives the van der Waals coefficient Cg — 253 kRy aQ.E3 
Earlier, simpler methods have been used to estimate the 
polarizability and the van der Waals coefficient. A^um- 
mation of C-C interactions gives Cg = 200 kRy aQ,c3 and 
for a calculation of dipole modes usiac a discrete dipole 
model the resuh is Cg 350 kRy a^Musing a(0) = 570 
ap (experimental value from Ref. ^) we get Cg — 302 
kRy aQ, a result which lies in the same range as those 
from the other calculations. 

In Tables || and || the characteristic frequencies, mq, 
corresponding to London's empirical formulate 

where A and B denote the two fragments, are also given. 
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This formula provides an easy way of estimating isotropic 
van der Waals. coefficients for mixed pairs of atoms and 
molecules BZl'Ea 

In Table ^ we have only given the isotropic dispersion 
coefficients for the molecules, that is, we have used the 
averaged polarizability a(m) in Eq. (g), where 

a{iu) = ^[a^r,(iu) + ayy{iu) + a^^(m)] . (12) 

It is easy, however, to calculate the anisotropic correc- 
tions, in addition. A simple example is for two inter- 
acting identical linear molecules, where the anisotropic 
coefficients Cg and Cg control the orientation-dependent 
part of the long-range interaction according toEj 

Evdw{R, (^B, 4>A, 4>b) — 

'" An 

Ce + C'e {P2icos9A) + P2(cos0s)} + y ^'x 



(3-|m|)ir(^^A,0A)l"2""(0B,'^i3) 



R- 



where 9a {Ob) is the angle between the vector R from 
the center of molecule A to the center of B and the axis 
of molecule A (B). The other angle ((/>b) describes 
the rotation of molecule A (B) about R. With 



these coefficients can be writtenil 

-j^ poo 



CL = — I duai{iu)Aa2{iu) 



and 



^6 =^ 



du Aai(m) Aa2(ju) 



(14) 



(15) 



(16) 



Calculating the anisotropic coefficients for H2 we obtain 
C^/Cg = 0.08 and C^/Cq = 0.007. Accurate values are 
C^/Ce = 0.1 and C^/Ce = O.Ol.BThe anisotropy is thus 
underestimated slightly with the simple cutoff scheme de- 
scribed above. In Figures || and || our calculated azzi^u) 
and axxiiu) for H2 are compared with accurate results. 

For an atom or molecule a distance d outside apirface, 
the asymptotic van der Waals energy is given byE3 



vdW — 



id ~ zor ' 

where the van der Waals coefficient is 

r 1 -(■ Miu)-l 

An Jq eb[iu) + 1 

and with the van der Waals plane 

1 



(17) 



(18) 



Zo 



AnC- 



3 Jo 



du a(iu) \ ^ — d(iu) 

^ Ubitu) + 1 eb{iu) + 1 ^ ' 



In these expressions, eb(m) is the bulk dielectric function, 
and d(iu) is the centroid of the induced surface charge 
caused by an electric field oriented perpendicular to the 
surfac fl and varying in time like e"*. Our earlier calcu- 
lationsll3ll3 of C3 and Zq have used the exact electro- 
dynamics, Eq. (|5|), for the surface but the approximate 
treatment, Eq. (ph, for the atom and molecule. In this 
paper we use Eq. (^) also for the latter. In Table III 
calculated values for C3 and Zq are given for He and 
H2 outside jellium, showing a very good agreement with 
other more elaborate calculations. 

Including the orientational dependence that results 
from the anisotropy of the molecular polarizability, the 
energy for a honionuclear diatomic molecule is to first 
order given bytil 



(0) 



P2(C0S( 



(20) 



where Q is the angle between the molecule axis and the 



(13) surface normal. C\ is given by Eq. (p^) and 



du Aa{iu) 



ebjiu) - 1 
eh(m) + 1 



(21) 



In Table |l| the ratio V^s"^ 

is given for II2 outside 



jellium. Wc find this ratio to be around 0.05 in agreement 
with Ref. |l[ 

We have in this paper refined the electrodynamical 
treatment of atoms and molecules within our previously 
presented density functional framework, thereby unifying 
our approaches for objects of different sizes. The calcu- 
lated polarizabilities and van der Waals coefhcients are 
in good agreement with results in the literature. This 
gives a possibility to easily calculate these quantities for 
complex systems with useful accuracy. 
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TABLE I. Van der Waals coefficients Ce for pairs of identi- 
cal atoms (Ry atomic units). The static polarizabilities used 
for defining the cutoff are given in the second column (atomic 
units) and results from other calculations in the fifth one. 
The third column gives the frequency uq obtained from the 
London formula, Eq. dill). 
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Mo 


Ce 
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0.41 
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Mg-Mg 
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Ca-Ca 


154-'' 


0.25 


4430 


4010^' 




|37|, " Ref. 
a Ref. p. 



M '= Ref. 
"^Ref. Isgl 



Ref. 



Ref. 



TABLE II. Van der Waals coefficients Ce for pairs of iden- 
tical molecules (Ry atomic units). The static polarizabilities 
used for defining the cutoff are given in the second column 
(atomic units) and values from the literature in the fifth one. 
The third column gives the frequency mq obtained from the 
London formula, Eq. ([ll|). 
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38"= 
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200 fc'', 253 fc\ 350 fc'' 





Q(0) 




Ce 


H2-H2 


5.41" 


0.98 


21.5 


N2-N2 


11.77'' 


1.44 


149 


CO-CO 


13.T= 


1.37 


176 


HF-HF 


5.52"^ 


1.77 


40.4 


H2O-H2O 


9.64^ 


1.40 


97.4 


Ceo-C'eo 


5703 


1.24 


302 fc 




|29|, " Ref. |0|, '= Ref. |41j, Ref. |2|, 

^ Mean value of estimates from Ref. 
3l 



= Ref. |3L_^ 
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TABLE III. The van der Waals coefficient C3 and the van 
der Waals plane position Zq (Ry atomic units) for He and H2 
interacting with jeiiium. For H2 aiso the ratio between the 



anisotropic coefficient C- 



(2) 



and ' 



IS given. 
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0.59 


0.64* 



Ref. 



Ref. 



FIG. 1. Calculated van der Waals coefficients Cq (Ry 
atomic units) for all possible pairs of the atoms in Ta- 
ble |l| pjlattjc^-pgainst corresponding values from other calcu- 
lationsElt 
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FIG. 2. Our calculated a{iu) for He, compared with a more 
accurate calculation.tj 
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FIG. 4. Our calculated|-«22(m) for H2, compared with a 
more accurate calculation.Lj 
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